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We present observations of the dwarf galaxies Draco and Ursa Minor, the 
local group galaxies M32 and M33, and the globular cluster M15 conducted with 
the Whipple 10m gamma-ray telescope to search for the gamma-ray signature of 
self- annihilating weakly interacting massive particles (WIMPs) which may con- 
stitute astrophysical dark matter (DM). We review the motivations for selecting 
these sources based on their unique astrophysical environments and report the 
results of the data analysis which produced upper limits on excess rate of gamma 
rays for each source. We consider models for the DM distribution in each source 
based on the available observational constraints and discuss possible scenarios 
for the enhancement of the gamma-ray luminosity. Limits on the thermally av- 
eraged product of the total self-annihilation cross section and velocity of the 
WIMP, (cf), are derived using conservative estimates for the magnitude of the 
astrophysical contribution to the gamma-ray flux. Although these limits do not 
constrain predictions from the currently favored theoretical models of supersym- 
metry (SUSY), future observations with VERITAS will probe a larger region of 
the WIMP parameter phase space, {av) and WIMP particle mass (m^). 

Subject headings: gamma rays: observations — dark matter 



1. Introduction 



The existence of dark matter (DM) is supported by a v ariety of observatio nal data in- 
cluding measurements of the cosmic microwave bac kground (ISpergel et al.l 12007). the large- 



scale distribution of galaxies (ITegmark et al.ll2004j ). and gravitational lensing (IClowe et al. 
20061 ). In the ACDM cosmological model that is currently favored by t hese data, DM cora - 
prises approximately ~26% of the total energy density of the universe (ISpergel et al.l 120071 ). 



However, the nature of the particles that constitute DM remains unknown. A popular DM 
candidate is weakly interacting massive particles (WIMPs), which existed in thermal equi- 
librium during the early universe and later decoupled as the universe expanded. Since the 
time of decoupling, the WIMPs have remained non-relativistic, behaving as a coUisionless 
fluid on all but perhaps the shortest spatial scales. In order to reproduce the observed rehc 
density of DM, this hypothetical particle would need to have a cross section on the scale of 
weak interactions. A stable particle with these properties, the lightest neutralino x-, can be 
accommodated in theories of supersymmetry (SUSY). 

The mass of the neutralino is const rained to be > 6 GeV by CMB measurements and 
accelerato r searches (Bottino et al. 20031) and < 100 TeV by the unitarity limit on the ther- 
mal relic (iGriest fc Kamionkowskil Il990l ). In the conventional SUSY scenarios, the neu- 
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tralino is a Majorana particle which can efficiently self-annihilate in astrophysical environ- 
ments with high DM density producing secondary particles including high-energy gamma 
rays. The former and current generation of Air-Cherenkov Telescopes (ACTs), including 
Whipple, HEGRA, CANGAROO-III, VERITAS, H.E.S.S., and MAGIC are sensitive in the 
gamma-ray energy range from below 100 GeV to above 10 TeV and can therefore make a 
substantial contribution to the search for the signatures of DM self-anni hilation. Recently 



several ACTs have detected gamma rays fr om the Galactic Center (G.C.) (lAlbert et al.ll2006 



Aharonian et al.ll2004l : iKosack et al.ll2004j). Although a more t raditional astrophysica l origi n 



for this signal is currently favored (lAtoyan &: Dermerl l2004l : lAharonian fc Neronovl l2005l ) . 
DM self-ami ihilation has been proposed as a possible explanation for these observations 
flHornsl[2005l ). 

We present observations taken with the Whipple 10m telescope of five astrophysical 
sources with the purpose of detecting the signature of DM self-annihilation. Section [2] sum- 
marizes the motivations for selecting each source. In Section [3] we discuss the signature of 
DM self-annihilation into gamma rays and its dependence on the source astrophysics and the 
particle physics properties of the WIMP. In Section H] we review the atmospheric Cherenkov 
technique and the methods used to analyze the data. Results of the data analysis are de- 
scribed in Section 14. 3[ Models for the DM distribution and scenarios for the enhancement 
of the gamma-ray flux are presented in Section [51 We conclude in Section [6] by discussing 
the implications of these observations for the parameter space of allowed SUSY models. 



2. Review of Observational Targets 

Because the gamma-ray signature of the neutralino is proportional to the square of the 
local density, the spatial scales that contribute to the total gamma-ray flux from a DM halo 
are much smaller than the spatial scales contributing to the halo mass. On these small spatial 
scales, the evolution of DM is typically driven by its interaction with baryonic matter, which 
dominates the gravitational potential. The influence of baryons in the form of dense stellar 
populations, molecular clouds, and central black holes could potentially lead to a much higher 
central DM concentration than that inferred from the large-scale DM distribution and thus 
substantially enhance the annihilation signal. Therefore it is attractive to consider sources 
that represent a diverse set of astrophysical environments, which could boost the gamma- 
ray luminosity. We have selected the dwarf galaxies Draco and Ursa Minor, the local group 
galaxies M32 and M33, and the globular cluster M15 for observations with the Whipple 10m 
telescope based on the analysis of observational data and various potential scenarios for DM 
enhancement in these objects. 
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2.1. Draco and Ursa Minor 



Dwarf spheroidal galaxies have attr acted considerable theoretical attention as poten- 
tial DM annihilation gamrr i a-ray sources (Balt z et al.l l2000l : iTylerl 120021 : IStrigari et al.l 120071 : 
Bergstrom &: Hooperl l2006l : IColafrancesco et al.. ,2007i) due to their large observed mass-to- 
light ratios (M/L). Studies of s urface brightness morphology have fo und a smooth symmetri- 



cal profile in the core of Draco (IPiatek et al.ll2002l: ISegall et al.ll2007f) but significant structure 



in the central regions of Ursa Minor jBellazzini et alibood ). which may be indicative of tidal 



interaction with the Milky Way. The observed stellar velocity dispersions in Draco and Ursa 
Minor imply that the dynamics of these systems are dominated by DM on all spatial scales 
and provide robust lower bounds on the astrophysical contribution to the gamma-ray flux 
in these systems. 

Both Draco and Ursa Minor possess low metall icity stellar populations with an age of 
~10 Gyr (lAparicio et al.l 1200 ll : IShetrone et al.ll200lh . The absence of recent star formation 
suggests that neither system has undergone a significant merger or accretion event since this 
early star formation epoch. Furthermore, because the two-body relaxation time in these 
galaxies significantly exceeds the Hubble time, it is unlikely that baryonic matter played 
a significant role in shaping the present-day DM distribution. Therefore primordial DM 
fluctuations formed on small spatial scales during initial violent relaxation could have been 
preserved and may boost the gamma-ray flux from these objects. 



2.2. M15 



Although there is no observational evidence for the presence of significant DM in globular 
clusters, the association of globular clusters and DM halos fits naturally into the standard 
pa radigm o f hierarchical structure formation. In the primordial formation scenario proposed 
by iPeebled ( 1l984j ). globular clusters are formed in DM overdensities in the early universe 
and may therefore retain a significant fraction of this primordial halo in the current epoch. 
Given that the extremely dense stellar cores of globular clusters dominate the gravitational 
potent i al of t hese systems, the observable effects of an extended DM halo may be minimal. 
Moord (Il996l ) argued that the presence of tidal tails in some globular c lusters suggests that 
globul ar clusters are not em bedded in DM halos. Recent simulations (IMashchenko fc Sills 
2005al lbl: | Saitoh et al.l l2006l ) have challenged this picture, showing that an extended halo 
may be compatible with the observable properties of globular clusters, although much of the 
original halo mass could be stripped by tidal interactions with the host galaxy. 



The proximity and potentially high central DM density of M15 favor this source for 
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indirect DM s earches. With a core radius of ~0.2 pc and extreme central density in excess of 
10^ M0 pc~^ (IDuU et al.lll997l ). M15 is the prototype for the core-collapsed globular cluster. 
During core collapse, the globular cluster is predicted to relax through stellar two-body 
collisions to a power-law d ensity profile that extends down to the smallest observable scales 
(IBinney fc Tremaindll987l ). If M15 was originally embedded in a DM halo, this evolutionary 
process must significantly compress the central DM distribution and dramatically enhance 
the gamma-ray flux. However, the poorly understood process of kinetic heating of DM in 
the core of the cluster by stars and hard binaries could lead to a depletion of DM from this 
region. 



2.3. M32 



The model for compression of DM through the gravitational contraction of dense stellar 
populations may also apply on the galactic scale. M32 is the closest compact elliptical 
galaxy and may have for med in a merging event between M31 and a low-luminosity spiral 
galaxy (IBekki et al.ll200ll ) in which the disk component of M32 was tidally stripped. Stellar 
kinematical data strongly support the presence of a si ngle supermassive compact object in 
the center of the galaxy with a mass of 2-4 x 10^ M© (j Joseph et al.ll200ll ). The core of M32 
has a r elatively homogene o us stellar population w i th an interna e diate age of approximately 
4 Gyr JCorbin et al.lbnOll : Idel Burgo etaPbonil ). Ihauer et all JlQQsl ) estimate M32's core 
relaxation time scale to be 2-3 Gyr, implying that the nucleus of M32 had at least a few 
relaxation times to evolve since the last significant merging event. Such events in which a 
massive black hole binary is formed are predicted to deplete the ce ntral density by evacuatiri g 
stars and destroying any potential DM cusp in the galaxy core ( iMilosavljevic et al.ll2002l ). 
Collisional two-body relaxation of a stellar population around a black hole is analytically 
predicted to result in a steady state power-l aw stellar density profile with power-law index 
between 3/2 and 7/4 (IBahcall &: Wolflll976l ). Optical and infrared data indicate a stellar 
density pro file compatible with a power-law index in the range 1.4-1.9 at the resolution limit 
of 0.07 pc (ICorbin et al.ll200ll : iLauer et al.lll998l ). Because the condensation of baryons in 
galactic nuclei may greatly enhance the central concentration of DM halos, the stellar density 
i n the core of M32 , in excess of 10^ Mq pc~^ and the highest known among nearby systems 



(iLauer et al.lll998r ). makes it a promising candidate for the detection of DM annihilation. 
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2.4. M33 



By observing astrophysical systems capable of rapid evolution, one may be able to 
overcome dynamical limitations on the neutralino annihilation rate, if it is limited by the 
scattering of WIMPs into a very small annihilation region in the galactic nucleus. M33 is 
remarkable for the small relaxation time, ~3 Myr, in its stellar nucleus of approximately 
0.2 pc, which results from the high stellar d ensity , 5x10^ pc~^, and extremely low 
velocity dispersion, 21 km s~^, in this region (ILauer et al.l 119981 ). M33 is a low-luminosity, 
DM dominated, bu 



Corbelli fc Salucci 



( iGebhardt et al. 



jeless spiral galaxy with a dark halo mass of approximately 5.1x 10 M© 
20001 ). The mass of the black hole in its center is less than 1.5x10^ Mq 



2001 



Merritt et al.ll200ll ). The stellar population in the nucleus of M33 can 



be modeled by two bursts of star formation 2 and 0.5 Gyr ago suggesting the possibility of 
a merger in the last ~Gyr. However, due to its rapid collisional relaxation time, M33 could 
have developed a core-collapsed nucleus in the period since the last merging event. 



DM Annihilation Flux 



The differential flux of gamma rays from WIMP annihilation along a line of sight is 
given by 

d(j){ip,AQ) _ (av) /dN{E,m^ 



dE 



dE 



d^l / (s,-ip,^l] ds, 



V. 



where p is the DM mass density, is the mass of the WIMP particle, (av) is the thermally 
averaged product of the total self-annihilation cross section and the velocity of the WIMP, 
dN{E,m^)/dE is the differential yield per annihilation, AQ is the solid angle observed, and 
ip is the direction of the line of sight integration. 



This expression can be factored as 

d(j){ip,An) _ 
dE ~ 



{av) 



3 X 10-26 cm3 s-i 



100 GeVV fdN{E,m^)/dE\ J{^f,An) 



■)( 



(2) 



10-2 GeV^ 



1.45 X 10^' 



where 



^1% 



6.64 X 10-^2 cm"2 g-i jg i% Qf ^^^q integral Crab Nebula flux above 100 GeV as 
extrapolated from the power-l aw fit of 3.2 x 10~ ^ ^ (E/ TeV)"^'^^ cm-2 s~^ TeV~^ reported by 
Hillas et al.l (119981 ). Following iBergstrom et al.l (Il998l ). the astrophysical component, which 
depends on the DM density profile, is expressed in terms of the dimensionless quantity J, 

1 



J{^p,An) 



pIRh 



An 



dQ / p [s,ip,Q \ ds 



(3) 
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which we normalized to the critical density pc = 9.74 x 10 g cm ^ and the Hubble radius 
Rh = 4.16 Gpc. 



3.1. DM Halo Profiles 



Estimation of J for a particular astrophysical object critically depends on the DM 
density profile. Numerical cold dark matter (CDM) simulations, applicable in the regions 
where DM dominates the overall gravitational potential, indicate the existence of a universal 
density profile across the spectrum of halo masses, from dwarf galaxies to galaxy clusters, 
which can be fit by. 



p(r) 



r 

Ps \ 



■ I 



(4) 



Simulations have consistently found that the large-scale (r ^ r^) asymptotic behavior is 
compatible with /5~3. The value of the inner logarithmic slope, 7, is less certain due to 
numeric al resolution effects at the smallest scales. The so-called NFW profile, (a,/3,7) = 
(1,3,1), (INavarro et al.lll996l ) has an inner asymptotic of r~^. P rofiles with a steep er inner 
power-law cusp, (a,/?,7) = (1.5,3,1.5), have also been suggested (IMoore et al.lll999l ). Recent 
high-resoluti on simulations have p ointed to an intermediate value for the inner power-law 



cusp 7 '-^1.2 (jDiemand et a. 



scales ( INavarro et al 



2004 



2005n and a possib le asymptotic shallowing approaching smaller 



Merritt et al.ll2006f ). 



The density profiles of simulated DM halos can described by their virial mass niyir and 
concentration c, which ar e related to p.^ and r» . Using a sample of simulated DM halos with 
masses 10^^-10^^ M0. lBullock et al.l (120011 ) found that the median halo concentration at 
2; = is correlated with virial mass and can be well approximated by the expression. 



9 



1.5 X lOi3/i-iM0 



-0.13 



(5) 



with a scatter of Alogc ~ 0.14. 



Observations of low surface brightness galaxies (Ide Blok et al.ll200ll ) have motivated an 
alternative form for the halo density profile which fits the rotation curves of these galaxies. 



Pir) = Ps 



r 

1 + - 



1+ I - 



_/3-7 



(6) 



A specific choice of (a,/?,7) = (2,3,1) is known as the Burkert profile (jBurkertlll995l ). Given 
the present state of uncertainty regarding the inner shape of DM halos, we use both the NFW 
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and Burkert profiles as representative of tlie possible range of inner density asymptotes in 
DM dominated fialos. 

The existence of DM halo substructures is a central prediction of CDM cosmology and 
may significantly enhance the DM annihilation flux as compared to that predicted for a 
smooth halo. Numerical simulations have predi cted that these substr uctures have a power- 
law mass function dN/dM ~ M~" with a ~ 2 (iDiemand et al.ll2006l ). The mass spectrum 
cutoff, mo, for these structures is set by the free streaming and coUisional damping length 
scales in the early universe and depends on the microphysical properties of the neutralino. A 



neutralino with mass = 100 GeV is estimated to have mo ~ 10 ^ Mq (iGreen et al.ll2005l ). 



The abundance of these substructures at 2; = is uncertain and depends on the fraction 
that survive disruption during the hierarchical merger and accretion processes. Current 
CDM simulations do not have sufficient resolution to model explicitly the mass function 
at all relevant scales. A recent estimate of the contribution of substructure to the DM 
annihil ation flux places a lo w er limit on the enha ncement to J relative to a smooth halo 
of 2-3 (IDiemand et al.l 120071 ). IStrigari et al.l (120071 ) obtain an upper bound of ~100 on the 
substructure enhancement by summing the contribution of all substructures above mo — 
10~^ Mq and assuming a mass function with tt=1.9. 

In the galaxies in which baryonic matter dominates the mass profile in the central 
region, the observational constraints on the astrophysical enhancement J are weak. Due to 
the large concentration of baryonic matter on small scales, the evolution of DM in galactic 
nuclei is likely significantly affected by the gravitational interaction of neutralinos with a 
central black hole, high density stellar populations, and gas. Evolutionary condensation of 
baryons in the core of a DM halo gives rise to a DM density enhancement beyond what 
would be expected from the gravitational interaction of DM alone. Using CDM simulations 
that combine d both dissipa t ionles s DM particles and a baryonic component composed of 
gas and stars, iGnedin et al.l (120041 ) observed that baryon condensation increases the central 
concentration of a DM halo. 

The adiabatic compression of DM in the core of a halo through the slow growth of 
a central black h ole has been suggeste d as a mechanism for increasing the flux from DM 
self-annihilation (iGondolo fc Silkl Il999l ). However, the magnitude of this effect depends 
strongly on the ratio of the initial and final masses of the black hole, its initia l alignment 



with t he center of the DM ha. lo, and the merging history of the galactic nucleus (jUllio et al. 



2OOII ). iMerritt fc Cruzl (120011 ) have shown that a typical merger event between black holes 
of comparable mass destroys cuspy profiles. The effect of gravitational scattering of DM 
particles by infalling baryons in the central regions of galactic nuclei may also be significant. 
The transfer of momentum to the dark matter component in these interactions should lead 
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to the partial evaporation of DM from the center of the galaxy. The ejection of DM particles 
through the gravitational slingshot process in the vicinity of binaries will further enhance 
the DM outflow from the centers of galactic nuclei. It might be expected that for a variety of 
astrophysical objects the main contribution to J originates from the regions of high baryonic 
density. The diversity of conditions affecting the DM distribution in the core of a galactic 
nucleus (presence or absence of a black hole, stellar velocity dispersion, density of baryonic 
matter, merging history, etc.) argues for observations of a variety of astrophysical objects 
in which DM self-annihilation might be amplified due to the gravitational interaction with 
baryons. 



3.2. DM Annihilation Spectrum 



The differential yield of gamma-ray photons per neutralino self-annihilation is a sum 
over final-state contributions, 



dN{E,m^ 
dE 



2b. 



+ 



E 



9 



dNi{E,m^ 
dE 



(7) 



where bx indicates the branching fraction of neutralino self-annihilation into a specific final- 
state channel X. The first two terms represent annihilations into mono-energetic pho- 
tons through either one or two photon channels. The last term is a sum over all chan- 
nels which contribute to the continuum fiux, which arises primarily from the decay of vr*^ 
me sons produced in the hadronization of the fermion and boson final states (for a review 



see 



Jungman et al.l (119961 )). The two-body annihilation modes into mono-energetic photons. 



77 and Zj, are significantly suppressed as compared to the continuum com ponent, with 
branching ratios ~10^^-10~^ (iBergstrom &: UUiol 119971 : iBergstrom et al.lll998l ). 



The differential spectrum of the 7r° decay component is relatively featureless and sim- 
ilar for all channels. It falls exponentially at high energies terminati ng at where it i s 
enhanced by internal bremsstrahlung from charged decay components (IBirkedal et al.ll2005l ). 



resulting in an edge-like feature at = m^ (see Figure [T]). Decays into all quark and 
bosonic states differ only slightly in the amplitudes of the vr'' and internal bremsstrahlung 
components. However, decay into r leptons generates a significantly harder spectrum, due 
to direct production of 7r° mesons in processes such as — > vr^vr^z/. In this work, the bb 
and r+r~ spectra are used to contrast the detection prospects in the case of a soft or hard 
neutralino self-annihilation spectrum, respectively. 



Although observing regions of high baryonic density improves the chances of indirect 




Fig. 1. — Spectral energy density per annihilation for a neutralino of mass 500 GeV anni- 
hilating to bb (solid line) and r"^r~ (dashed line). These spectra were generated using the 
PYTHIA Monte Carlo package Jsiostrand et al.lboOlj ). 



neutralino detection, the high density of baryonic matter comes with the price of potentially 
high astrophysical gamma-ray backgrounds. High-energy processes that take place in galac- 
tic nuclei, such as acceleration of particles in supernova shocks or the jets formed by accreting 
black holes, interactions of cosmic rays with molecular clouds, etc., can generate detectable 
gamma-ray fluxes from these regions. In these processes, particles producing gamma rays are 
most likely accelerated stochastically resulting in a power-law differential energy spectrum 
with index ~2. The truncation of the spectrum of gamma rays from neutralino annihilation 
at the neutralino mass can, in principle, be mimicked by traditional astrophysical processes. 
For example, local absorption of gamma rays through pair production in high-density op- 
tical and infrared diffuse photon fields would cause an exponential cutoff in the spectrum. 
However, such truncation would depend on the specific properties of the source and would 
not be a common observational feature across multiple objects. 



4. Observations and Analysis 

The TeV gamma-ray observations of the five sources reported in this work were taken 
with the Whipple 10m telescope located at the Fred Lawrence Whipple Observatory in south- 
ern Arizona, USA. The telescope is equipped with a camera consisting of 379 photomultiplier 
tubes, covering a field of view of diameter 2.4°, which detects the short-duration Cherenkov 
light flashes emitted by secondary particles generated in cosmic- and gamma-ray-induced at- 
mospheric cascades. A detailed description of the telescope optics and camera configuration 



is presented elsewhere (IKildea et al.ll2007l ). 



4.1. Observations 

The five sources were observed over the course of four observing seasons using two 
different modes of operation: ON/OFF and TRACK. The ON/OFF mode is characterized 
by a sequential pair of 28 minute runs in which the ON run is obtained by pointing the 
telescope in the direction of the source and the OFF run is taken at the same azimuth 
and elevation but offset by 28 minutes in right ascension for background estimation. The 
ON/OFF technique minimizes systematic errors due to the changing state of the atmosphere 
and variations in the night-sky brightness. The TRACK mode observations are not followed 
by dedicated OFF observations. Instead, a contemporaneous but unrelated OFF run of 
similar elevation and average night-sky background noise is selected and then analyzed with 
the TRACK mode observation in the same way as the ON/OFF pair. To reduce systematic 
errors due to variations in the night-sky background between ON and OFF fields, artificial 
noise was inje cted on a pixel by pixel basis into either ON or OFF runs to equalize them 



(IWeekeslll996l ). The total accumulated exposures in ON and TRACK modes and observation 



epochs for each source are summarized in Table [TJ 



4.2. Data Analysis 



The light dis tributions of the Cherenkov images were parameterized using a standard 
moment analysis (jHillaslll985l ). The parameter size characterizes the total amount of light 
in the image while the parameters width and length, distance, and a describe the shape, 
location, and or ientation of the imag e, respectively. A set of selection criteria, referred 
to as Supercuts ( iReynolds et al.lll993l ). were used to select gamma-ray-like images with an 
efficiency of ~50%, while rejecting >99% of cosmic rays. 



For events passing the Supercuts criteria, a histogram of the a parameter is generated. 
Gamma-ray events originating from a source in the center of the field appear as an excess in 
the region of the a-histogram with a < 15°. The residual background in the signal region of 
the ON-source observation is estimated from the number of events obtained from the same 
region of the OFF-source observation rescaled by the ratio of the OFF to ON livetimes. In the 
case of the TRACK mode observations, the background estimate from the matching OFF- 
source observation is scaled by an additional factor, equalizing backgrounds in a sideband 
region 20° < a < 65° on a run-by-run basis. For each source, the significance of an excess 
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in the signal region is evaluated with the maximum likelihood method of iLi fc Mai (Il983l ). 



The method of analysis described estimates the integral gamma-ray excess which is 
the convolution of the effective area of the Whipple 10m telescope shown in Figure [2] with 
the spectrum of the source. The gamma-ray collecting area is nearly constant at ~4 x 10^ 
cm^ for photon energies above 1 TeV and declines rapidly in the energy regime below ~400 
GeV. The differential gamma-ray detection rate, the product of effective collecting area and 
source spectrum, is shown for a Crab Nebula-like spectrum in Figure [31 The peak of the 
detection rate is ~400 GeV and is a relatively weak function of the index of the power- 
law spectrum in the range ~2-3. Despite the fact that the continuum neutralino spectrum 
can not be approximated as a power law, for neutralino masses exceeding approximately 
400 GeV the differential gamma-ray detection rate peaks in the vicinity of this energy as 
illustrated in Figure [3] for = 1 TeV. Thus the integral constraints discussed in this paper 
are appropriate for a source with a neutralino-like spectrum as they have been optimized to 
produce the best sensitivity for the sources which peak in the energy range 300-400 GeV. 
The integral method remains applicable for neutralino masses below 400 GeV. However, the 
sensitivity deteriorates rapidly as the peak of the detection rate falls significantly below 400 
GeV. 

In this analysis we have not searched for the monoenergetic line feature of the neutralino 
self-annihilation spectrum. The Whipple 10m telescope has an energy resolution of ~30% 
which ultimately limits the sensitivity to monoenergetic photons. The search for the contin- 
uum component of the spectrum can provide an equivalent or better sensitivity to neutralino 
annihilations as compared to searching for the monoenergetic line given that under the pre- 
ferred particle physics scenarios the branching ratios to two and one photon channels is less 
than one percent. A dedicated search for monoenergetic photons could be a goal of a refined 
analysis if an integral excess is detected. 

The accuracy of the reconstruction of the arrival direction of individual photons for the 
Whipple 10m telescope for a source in the center of the field of view is a function of photon 
energy and changes from ~20' at 300 GeV to ~2' at 10 TeV. The event selection criterion 
a <15° corresponds to an angular cut in the image plane of ~15' which translates into a 
solid angle AQ of 6 x 10~^ sr which is used in the calculation of J (Equation [3]). 

The Crab Nebula is the standard calibration source for ground-based gamma-ray as- 
tronomy. Observations of the Crab Nebula were taken during the same four year period as 
the observations of the sources of interest. These data were analyzed with the same tech- 
nique to calibrate the gamma-ray detection efficiency of the Whipple 10m telescope. The 
differentia l flux of the Crab nebula at 400 GeV is estimated using the spectral parameter- 



ization of iHillas et al.l (119981 ). By scaling the observed gamma-ray rate from the putative 
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DM sources to the contemporaneous rate of the Crab Nebula, the systematics due to the 
changing performance of the telescope optics and camera system were corrected for. 



4.3. Results 



Using the Supercuts criteria, the excess gamma-ray rate and its significance for each 
source was calculated under the hypothesis of a point source in the center of the field of 
view. No significant excess was detected from any of the five sources observed (see Table 
12]). In the absence of a significant detection, the 95% C.L. upp er limit on the rate from 
each source was calculated following the method of iHelend (Il983l ). The upper limit on the 
differential spectral energy density E'^dF/dE at 400 GeV is calculated for each source by 
scaling the gamma-ray rate to the observed Crab rate during the same observation epoch. 

The 95% C.L. upper limits for all sources are 6-9% of the Crab Nebula flux with the 
exception of M15 for which the exposure time was sig nificantly shorter. ACT observations 
of M32 were reported by the HEGRA collaboration (lAharonian et al.ll2003l ). This work 
places a 99% C.L integral flux upper limit of 4.4% of the Crab Nebula rate for M32 which 
is compatible with the results presented in this work. 



5. Analysis of Astrophysical Enhancement Factors 

The primary difficulty in constraining the parameter space of allowable SUSY models 
is due to the significant uncertainty in the astrophysical enhancement factor, J (Equation 
|3]). A lower bound on J can be estimated by extrapolating the DM density measured on 
large spatial scales (r ^ r^) into the small scales (where most of the DM annihilation signal 
originates) by using a profile with a weak cusp (NEW) or central core (Burkert). However, 
for each of the sources considered there may exist astrophysical mechanisms which could 
enhance the density of DM in the core of the halo and boost the luminosity due to neutralino 
self-annihilation by several orders of magnitude. A discussion follows for each source which 
presents both a conservative estimate of J as well as possible scenarios for its enhancement 
which take into account the source's unique astrophysical environment. 



5.1. Draco and Ursa Minor 



Because the gravitational potentials of Draco and Ursa Minor are dominated by DM 
on all observationally resolved scales, studies of stellar kinematics in these systems provide 
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Table 1. Summary of observation period and exposure time for ON and TRACK 

observing modes 



Source 


Period 


ON Exposure 
(hours) 


TRACK Exposure 
(hours) 


Total Exposure 
(hours) 


Draco 


Mar 2003 - Jul 2003 


7.4 


6.9 


14.3 


Ursa Minor 


Jan 2003 - Jul 2003 


7.9 


9.3 


17.2 


M32 


Sep 2004 - Dec 2004 


6.9 





6.9 


M33 


Oct 2002 - Dec 2004 


7.9 


9.2 


17.0 


M15 


Jun 2002 - Jul 2002 


0.2 


1.0 


1.2 



Table 2. Detected Gamma-ray rates, inferred upper limits on the gamma-ray rate, 
significances, and differential gamma-ray flux upper limits 





Excess 


95% C.L. Upper Limit 




at 400 GeV Upper Limit 


Source 


(7 min~^) 


(7 min~^) 


a 


(erg cm ^ s ^) 


(% Crab) 


Draco 


0.001 ± 0.066 


0.14 


0.02 


5.10 X 10-^2 


6.23 


Ursa Minor 


0.075 ± 0.070 


0.20 


1.07 


7.30 X 10-^2 


8.94 


M32 


-0.240 ± 0.170 


0.21 


-1.44 


6.00 X 10-^2 


7.34 


M33 


-0.012 ± 0.085 


0.16 


-0.14 


5.75 X 10-12 


7.04 


M15 


0.496 ± 0.269 


0.94 


1.80 


2.68 X 10-11 


32.8 
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robust constraints on their DM density profiles at radii greater than ~0.5 kpc. IStrigari et al. 



( 20071) have used t he radial velocity data sets compiled by IWilkinson et al.l (120041 ) and 
Munoz et al.l (120051 ) to derive constraints on the parameters and ps of their DM halos 
under the assumption that the DM follows an NFW density profile. The best-fit contours 
in the Tg-ps plane are shown in Figures HI and [6] for Draco and Urs a Minor, respectively. A 
similar analysis of the Draco data set by iMashchenko et al.l (120061 ) which considered both 
NFW and Burkert DM density profiles obtained similar constraints on and p^, shown 
in Figures HI and El Both analyses find a region of degeneracy in the r^-p^ plane which is 
attributable to the weak constraints on the stellar velocity anisotropy. The region of degen- 
eracy for these models, however, is nearly parallel to the isocontours of J which results in a 
much smaller uncertainty for this pa rameter. Conservat i ve all owab le ranges for J calcul ated 
using the mass models presented by IMashchenko et al.l (120061 ) and IStrigari et al.l (120071 ) are 
summarized in Table |31 

The inner logarithmic slope (7) of the DM density profile in Draco and Ursa Minor 
is not observationally constrained. A value of 7 > 1 could potentially enhance by 2-3 
orders of magnitude the value of J with respect to the estimates presented in Table |31 The 
presen ce of DM substructu res is another potential factor contributing to the enhancement 
of J. IStrigari et al.l (120071 ) calculated an upper bound of ~100 on the enhancement due 
to substructures for a generic DM halo. The distortion of the gravitational potential by 
DM substructures may be reflected in the distribution of stellar populations. Stellar lumps 
within the central 10 arc minutes of Ursa Minor have been detected, and the study of the 
stellar pr oper motion suggests that th e lifetime of these structures should be no longer than 
~5 Myr (lEskridge fc Schweitzeill200ll ). If proven statistically significant, these observations 
may indicate the existence of small scale substructures in the DM distribution. However a 
conventional astrophysic al explanation such as the projection of a cold extratidal population 
of stars is also possible (IWilkinson et al.l 12004 ). 



5.2. M15 



Models for the M/L of M15 are consistent with a purely baryonic mass profile (Ivan den Bosch et al 



20061 ) and therefore indicate that the mass of a putative DM halo must be significantly less 
than the mass of the stellar component within the observable extent of the cluster. How- 
ever, given the compact nature and extreme central density of M15, this constraint is not a 
significant limitation on the potential DM annihilation signal, as evident from the following 
estimates. 

Consider a hypothetical halo model for M15 parameterized by a virial mass nimr and 



- 16 - 



concentration parameter c. The virial mass is unlikely to be significantly less than estimates 
of the present baryonic mass of the cluster of ~5x 10^ Mq ( IDuU et al.lll997l : iMcNamara et al. 



2004 ). Conversely, because of its distance from the Milky Way, the halo mass can not be 



significantly greater than ~10^ Mq, or dynamical friction would have resulted in the inspiral 
of the cluster in less than a Hubble time. Based on these constraints, we adopt a range for 
rrim r of 5x10^ Mp-5xl0 '' Mq and consequently a range for c of 32-82 as estimated with 
the iBuUock et al.l (120011 ) relationship. Equation [5 Although the c (m„ir) correlation was 
derived with simulated halos of mass ~10^^-10^^, IColm et al.l (120041 ) have found that this 
relationship holds for simulated halos with masses down to ~10'' Mq. To account for tidal 
disruption due to the interaction with the Milky Way, we truncate the hypothetical model 
of the DM profile at the optical extent of the cluster, ~30 pc. 

For the adopted ranges of rrimr and c, we estimate J to be 7-150 if the DM halo 
follows an NFW profile. However, the adiabatic compression of DM in the core of M15 
co uld further enhance the anni hilation signal by several orders of magnitude. Simulations 
by iMashchenko fc Sillj (j2005bl ) of a two-component globular cluster with stars and DM 
have demonstrated an enhancement to the central DM density as the baryonic mass pro- 
file becomes more centrally concentrated through two-body interactions. To evaluate the 
enhancement to the central DM d ensity for our hypothetic al model of M15, we have used 
the adiabatic compression model of iBlumenthal et al.l (119861 ). With the assumption that the 
DM travels on circular orbits, this model relates the initial and final baryon and DM mass 
profiles by the equation, 

[MDMjirf) + Mhjirf)]rf. ^ ' 

The initial distribution of DM is described by an NFW profile parameterized by rrimr and 
c, while the initial baryonic mass profile is assumed to follow the DM mass profile with the 
cosmological baryon-to-DM ratio of 0.2, following the assumption that globular clusters are 
among the oldest gravitationally bound systems with cores that have not been significantly 
infiuenced by merger events during their evolution. For the final baryonic mass distribution, 
we adopt a cored profile with (a,/5,7) = (2,2.6,0), = 0.04 pc, and p,^ ^10^ M^r^ pc ~'^ which 



approximates the nonparametric stellar mass profile presented in iGebhardt et al.l (119971 ). 
The predicted profile after adiabatic compression, hereafter denoted as NFW+AC, is shown 
in Figure [7]for the case of rrimr = 10^ Mq and c = 50. The density of DM interior to ~10 pc 
is enhanced by a factor ~10-10^ resulting in an increase in J of ^10^. The truncation of the 
DM halo at 30 pc has a negligible effect on the total annihilation signal. Figure [S] illustrates 
the scaling of J with the assumed virial mass and concentration of the DM halo. The range 
of J for the models with and without adiabatic compression are summarized in Table [31 
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5.3. M32 



The stellar surface brightness profile of M32 is characterized by a bulge with a half- 
light radius of 100 pc and a faint surface brightness e xcess beyond 3 00 pc which could be 
interpreted as the rem nant of a tidally stripped disk (iGrahaml |2002| ). A V-band M/L of 
2.51 interior to ~50 pc fivan der Marel et al.lll998l ) is consistent with the expectations of an 
intermediate age stellar population. Due to the absence of a significant disk component, the 
available kinematical data do not constrain the presence of an extended DM halo in this 
system. Its proximity to M31 has lead to speculation that M32 may be the remnant of a 
normal elliptical or late-type spira l galaxy which w as tidally stripped as it passed through 
the disk of M31 (IBekki et al.ll200ll : IChoi et al.ll2002l ). The extent to which tidal interactions 
may have stripped the DM halo in this event is unclear. 

The mass of the M32 DM halo can be estimated under the hypothesis that its progenitor 
was a late-type spiral by using the correla tion betw e en bu lge velocity dispersion Oc and 
maximur n circular velocit y Vc presented by iFerraresd (|2002| ) for a sample of nearby spiral 
galaxies. IFerraresd (120021 ) notes that may be related to the virial mass of the DM halo. 



using the correlation between them observed for simulated halos by iBuUock et al. 



(120011 ). By combining these two correlations the following relation between the bulge velocity 
dispersion and DM halo mass is obtained: 



1012M 



5.6 







Or. 



200 km/s 



2.79 



(9) 



The average bulge velocity dispersion in M32 of 76 ± 10 km s~^ ( Ivan der Marel et al.lll998l ) 
suggests a virial halo mass of 2.5-5x10^^ Mq. By further assuming that the DM halo of 
M32 follows an NEW or Burkert profile, the scale radius and normalization of the halo are 
fixed by a choice of the concentration parameter c, predicted to be in the range 10-20 by 
Equation [51 Constraints on the parameters for both NEW and Burkert DM density profiles 
are plotted in EigureO The estimated range of J for these mass models is found to be 1-10 
(see Table [3]). The models with the range of m^ir and c considered here are compatible with 
the the observed M/L in the interior of M32. 

With its extreme central density of 10'' M© pc~'^ and two-body core relaxation time of 
2-3 Gyr, M32 could possess an enhanced DM cusp produced through the compression of its 
DM density profile by baryonic infall and the growth of the central black hole. However, it is 
likely that M32 has undergone multiple mergers during its evolutionary history which could 
have significantly altered the DM density profile. The degree to which these effects may 
deplete the central density of DM in this system is uncertain, and therefore the estimate of 
J obtained with the adiabatic compression model should be considered an upper bound. We 
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constructed an adiabatic compression model using the same method apphed in the case of 
M15. The initial and final DM and baryon mass profiles for the case of m^j^ = 4 x 10^^ Mq 
and c = 13.7 are shown in Figure [TUl The range for J for the NFW and NFW+AC models 
is presented in Table El 



5.4. M33 



M33 is a late-type Sc galaxy with a substantial fraction of its baryonic mass in the form 
of neutral hydrogen. Me asurements of its rotation curve imply a lo wer limit for the mass of its 
dark halo of 5 x 10^° Mq JCorbelh fc Saluccilboool ) . ICorbellil J2003h have modeled the rotation 
curve, as derived from high-resolution velocity maps of CO, using a three-component density 
profile which includes a stellar nucleus of ~8xl0^ M0 interior to ~0.5 kpc, an exponential 
disk of mass ~3xl0^ M©, and an extended DM halo with a scale radius ~20 kpc. The data 
is inconsistent with a DM pro f ile wi th an intermediate power-law asymptotic as steep as r~^'^ 
as proposed by iMoore et al.l (119991 ) but can be well matched by either a Burkert or NFW 
profile. The best-fit region in the {ps, r^) plane for both NFW and Burkert p rofiles is shown 
in Figure [TTJ Ranges for J derived from the mass models of ICorbellil (120031 ) are presented 
in Table [31 

Observations of the core of M33 have shown evidence for a stellar nucleus with an 



effective radius of ~3 pc (iLauer et al.l Il99^ : [Stephens fc Frogelll2002l ) and a central velocity 
dispersion of ~24 km s~^ (jCebhardt et al. 2001). The exceptionally low velocity dispersion 
in this region sets an upper limit on the mass of a central black hole at ~1.5 xlO'^ Mq and 
is compatible with a relaxation timescale of ~10^ yr at 0.1 pc. Due to the unusually rapid 
evolutionary timescale in the stellar nucleus, M33 may be able to sustain a steep DM cusp in 
its core as DM could be rapidly replenished after a merger or accretion event. The kinematics 
of this region would be conducive to the growth of an intermediate mass black hole with a 
relatively large initial to final BH mass ratio and the creation of a DM cusp in the vicinity 
of the BH. The nucleus of M33 hosts the most luminous steady X-ray source in the local 
group which is a lso associated with a rad io source and similar to the galactic microquasar 
GRS 1915+105 JDubus fc Rutledgel l20o3 l . A small but significant time variability of 10% 
i n the X-ray lurainosit y and associated variability in spectral shape have also been observed 
Jla Parola et al.ll2003h . These data could be interpreting as supporting the existence of an 
accreting BH in the nucleus of M33. However, the exact enhancement to J through these 
processes is difficult to estimate given the unknown merger history and therefore an upper 
bound for J is not presented. 
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Table 3. A summary of estimates of J for each source 



Source 


Distance (kpc) 


Model 


J min 


J max 


Draco 


80 


Burkert 


1 


10 






NFW 


4 


40 


Ursa Minor 


66 


NFW 


4 


20 


M15 


10 


NFW 


7 


150 






NFW+AC 


8x10^ 


2x10^ 


M32 


776 


Burkert 


0.6 


9 






NFW 


1 


9 






NFW+AC 


2x10^ 


10*^ 


M33 


840 


Burkert 


0.03 


0.2 






NFW 


0.2 


0.6 



Note. — The contribution of DM substructure may poten- 
tially enhance the estimates of J shown here by a factor > 2 
and < 100 as discussed in Section [3?T1 
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6. Limits on SUSY Parameter Space 

For the prediction of gamma-ray fluxes from neutralino self-anniliilation, the framework 
of the minimal supersymmetric extension to the standard model (MSSM) was used. With 
several simplifying assumptions, the MSSM can be reduced to the seven parameters fi, M2, 
tan/?, mji, rriq, At, and Af,. Random sca ns of these parameters were performed utihzing 



the DarkSUSY code (iGondolo et ahl 120051 ) . and, for each model, the mass m^, (av), and 
relic density Vt^Mh? of the neutralino were calculated. MSSM m odels consistent wit h the 



±3 standard deviation bounds on VLouh? measured by WMAP (ISpergel et al.l 120071 ) were 
selected for comparison with the derived limits on (erf). Figure WI\ shows a projection in the 
{av)-m^ plane of MSSM models that satisfy both the WMAP constraint and the bounds 
placed by accelerator experiments. The majority of the models with neutralino mass above 
100 GeV are concentrated in the band with typical {av) ^ 3 x 10~^^ cm^ s~^ and extending 
to rriy^ ~ 2 TeV. 

For a typical choice of MSSM model parameters, the self-annihilation will predominantly 
proceed through some combination of the flnal states 66, tt, W~^W~ , or ZZ. The gamma-ray 
spectra of these channels are similar since they all result from the decay of neutral pions 
produced in the hadronization of the annihilation products. The r~*"r~ channel produces 
a signiflcantly harder spectrum as discussed in Section 13. 2[ Although a neutralino which 
annihilates predominantly to the r+r~ channel would improve the prospects for detection 
by an ACT, for the MSSM models considered here, the branching fraction of this channel is 
never more than ~10%. Three spectra generated using the PYTHIA Monte Carlo code were 
adopted to cover the possible range of dN{E,m^)/dE: the bb spectrum, the t~^t~ spectrum, 
and a composite spectrum with BR{xx ~^ bb) = 0.9 and BR{xx ~^ t~^t~) = 0.1. 

Table 4. Upper limits on (av) of the neutralino 







95% C.L. Upper Limit on 


{(jv) for = 


1 TeV (cm^ s-^) 


Source 


J 


bb 


r+r- 




Draco 


13 


< 1.9 X 10-21 


< 1.2 X 10- 


-22 


Ursa Minor 


9 


< 3.9 X 10-21 


< 2.5 X 10" 


-22 


M32 


3 


< 1.2 X 10-20 


< 8.0 X 10- 


-22 


M33 


0.4 


< 8.0 X 10-20 


< 5.2 X 10- 


-21 


M15 


32 


< 5.1 X 10-21 


< 3.3 X 10 


-22 
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Following Equation [21 the upper limit on {av) as a function of for a source with an 
astrophysical enhancement factor J and an upper limit on the detected rate of gamma rays 
(95% C.L.) is given by 



{av) 



3 X 10-26 cm3 s-i 

01% 



< R-y (95% C.L.) 



100 GeV 



2 1.45 X 10^ 
J 



^ ' V 10-2 Gev-1 ; 



(10) 



where the assumed form of neutralino annihilation spectrum is convolved with the energy- 
dependent effective area of the Whipple 10m telescope, A[E), shown in Figure El Table H] 
presents limits on {av) derived for a neutralino of mass 1 TeV annihilating through the hh 
and T^T~ channels. Figure [121 shows the limits on {av) as a function of for the NFW 
mass models of Draco and Ursa Minor and the adiabatic compression model (NFW+AC) of 
M15. Because the effective area of the Whipple 10m telescope rapidly dechnes below ~400 
GeV, the limits on {av) are most constraining for neutralino masses above this energy as 
discussed in Section 14. 2[ 

Using the most conservative estimates for J, the limits on {av) are 10^-10^ times greater 
than the range predicted for the MSSM models considered in this analysis. The DM mass 
models of Draco and Ursa Minor have the best observational constraints and the fewest 
uncertainties with regard to the unknown influence of baryonic matter and merging history. 
The lower limit on J for these galaxies is relatively insensitive to the assumption of a cusped 
versus cored DM density profile. The astrophysical contribution to the gamma-ray luminosity 
could be significantly enhanced if the effects of substructure or a density profile with an 
inner logarithmic slope > 1 were considered. Limits derived from ACT data on (av) and 
the mono energetic line cornpone nt {<Jv)^^ h ave previously been reported for observati ons of 
the G.C. dAharonian et al.lbooeh and M31 jAharonian et al.lboosl : iLavalle et allbood ). The 
limits reported by the H.E.S.S. collaboration derived from observations of the G.C. are among 
the most constraining with a 99% C.L. upper limit on (av) of 10-2^-10-2^. However these 
measurements come with a significant systematic uncertainty as they depend sensitively on 
the accurate modeling of the astrophysical background. For neutralino masses below 100 
GeV, EGRET data have also been shown to constrain (av) assuming a range of models for 
the distrib ution of DM in s ubstructures and the distribution of substructures in the Milky 
Way halo Jpieri et al.ll2007h . 



In order for the neutralino to be detectable by the Whipple 10m telescope, a significant 
enhancement of the DM density is required. Such enhancement may be consistent with the 
kinematics of M15, M32, and M33. Among the scenarios discussed, the adiabatic compression 
model for M15 provides the most quantitative estimate for J. In the potentially extreme 
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DM enhancement scenarios which may exist in M32 and M33, current hmits on the gamma- 
ray rate from these sources aheady have the potential to constrain the parameter space of 
allowed SUSY models. If one assumes that DM is composed of neutralinos with (av) ~ 3 x 
10~^^ cm^ s~^ and > 400 GeV, then the limits presented in this work rule out a value of 
J for these sources in excess of ~10^. 



7. Conclusions 

We have conducted a search for the gamma-ray signature of neutralino self-annihilation 
from five sources: the dwarf spheroidal galaxies Draco and Ursa Minor, the globular cluster 
M15, and the local group galaxies M32 and M33. Each of these sources was chosen as a 
favorable representative of different astrophysical conditions that could potentially enhance 
the neutralino density and gamma-ray self-annihilation flux. For a generic MSSM model of 
the neutralino, the self-annihilation flux is only detectable by the Whipple 10m telescope if 
such an enhancement is significant. 

A standard analysis of the data revealed no significant excesses, and upper limits on the 
gamma-ray flux from each source were derived relative to the flux of the Crab Nebula. We 
have derived limits on (av) of the neutralino as a function of its mass using models for the 
DM profile of each source and the generic differential gamma-ray spectrum of the neutralino 
self-annihilation constructed from two representative channels, bb and t'^t~. Using the upper 
limit on the gamma-ray rate measured from Draco and the most conservative estimate of the 
DM distribution in this source, we obtain 95 % C.L. upper limits on {av) of < 1.9 x 10~^^ 
cm^ s^^ and < 1.2 x 10~^^ cm^ s~^ for a neutralino of mass 1 TeV annihilating exclusively 
through the bb and t~^t~ channels, respectively. 

We have considered potential enhancements to the DM self-annihilation fiux including 
the effects of DM substructure and the adiabatic compression of the DM halo due to baryonic 
infall. These scenarios could enhance the annihilation fiux by as much as 10"^ and possibly 
higher. However the uncertainties of these estimates are large due to the poorly understood 
dynamics in the cores of the objects as well as their potentially complex merging histories. 
If one assumes that DM is composed of neutrahnos with {av) ~ 3x 10"^^ cm^ s~^, then 
some extreme enhancement scenarios for M32 and M33 may be ruled out. 

The current generation of ACTs such as VERITAS, H.E.S.S., MAGIC, and CANGAROO- 
III have the potential to improve significantly the sensitivity of these measurements and thus 
probe a larger region of the MSSM parameter space. For a source with a Crab Nebula-like 
spectrum, VERITAS has a fiux sensitivity ~10 times better than the Whipple 10m telescope 
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and a peak detection rate near 150 GeV. The lower energy threshold of VERITAS will al- 
low it to be sensitive to the low to intermediate neutralino mass range of 100 GeV-1 TeV. 
With the improved sensitivity of these instruments, conventional astrophysical backgrounds 
may become significant. The sensitivity of observations of the G.C to DM annihilations is 
already limited by the presence of such backgrounds. Observations of extragalactic sources 
such as those discussed in this work have the potential to avoid this limitation. Furthermore, 
the identification of the unique spectral signature of DM self-annihilations in two or more of 
these sources would effectively rule out a traditional astrophysical process. Next-generation 
ACT instruments such as the currently planned Cherenkov Telescope Array (CTA) and Ad- 
vanced Gamma Ray Imaging System (AGIS) will potentially be 10^-10^ times as sensitive 
as the Whipple 10m telescope and could perform dedicated deep observations with 10-10^ 
times longer exposure than the observations presented in this work. This may allow the 
exclusion of MSSM models with even the most conservative assumptions for the DM dis- 
tribution in the sources with the lowest anticipated astrophysical backgrounds such as the 
dwarf spheroidal galaxies. 
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Fig. 2. — Effective area of the Whipple 10m telescope as a function of energy after selection 
of gamma-ray-like events with the Supercuts criteria. The effective area below 150 GeV was 
assumed to be in all computations. 




Fig. 3. — Comparison of the differential detection rate of the Crab Nebula (solid hue) and 
a DM halo with (av) = 3 x lO^^e, rn^ = 1 TeV, and J = 10^ (dashed hue). The Supercuts 
selection criteria are used. 
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Fig. 4. — Constraints from stellar kinematics data on the parameters and ps of the Draco 
DM halo under the assump tion of an NFW profile. T he gray contour and solid c ircles indicate 
the best-fit mass models of IStrigari et al.l (120071 ) and iMashchenko et al.l (120061 ) . respectively. 
The parameters and ps are defined for the NFW and Burkert halo profiles in Equations H] 
and El 




Fig. 5. — Constraints from stellar kinematics data on the parameters and ps of the Draco 
DM halo under the assumption of a Bu rkert profile. The solid circles indicate the best-fit 
mass models of IMashchenko et al.l (120061 ). Thick solid lines indicate contours of constant J. 
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Fig. 6. — Constraints from stellar kinematics data on the parameters and ps of the Ursa 
Minor DM halo under the assumption of a n NFW profile. The grey contour indicates the 
best-fit mass models of IStrigari et al.l (120071 ). Thick solid lines indicate contours of constant 
J. 




r[pc] 

Fig. 7. — Shown is the M15 DM density profile before (thin solid line) and after (thick solid 
line) adiabatic compression modeled as described in the text. Thin and thick dashed lines 
show the same comparison for the baryonic density profile. 
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Fig. 8. — J for the NFW+AC model of M15 as a function of the concentration parameter 
of the initial DM halo. The solid and dashed lines show the scaling for a DM halo of total 
virial mass 5x 10^ M0 and 5x10^ M© respectively. 




Fig. 9. — Constraints derived using Equations [9] and [5] for M32 on the parameters and ps of 
its DM halo under the assumption of an NFW (left panel) and Burkert (right panel) profile. 
Dashed lines show contours of constant virial mass for 10^^ Mq (lower) and 10^^ M© (upper). 
Dotted lines denote the ±1 standard deviation bounds on the median halo concentration to 
virial mass relation, Equation [5l Thick solid lines indicate contours of constant J. 
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Fig. 10. — A comparison for M32 of the stellar density profile of Ivan der Marel et al.l (119981 ) 
(thick dashed line) with the modeled DM density profiles before (thin solid line) and after 
(thick solid line) adiabatic compression. The thin dashed line shows the assumed initial 
baryonic density profile. 




Fig. 11. — Constraints from CO rotation curve data for M33 on the parameters and ps 
of its DM halo under the assumption of an NFW (left panel) and Burkert (right panel) 
profile. The gr ay-shaded region indicates the 3 standard deviation constraints reported by 
Corbellil (120031 ). Dashed lines show contours of constant virial mass for 10^^ Mq (lower) and 



10^^ Mq (upper). Dotted lines denote the ±1 standard deviation bounds on the median 
halo concentration to virial mass relation, Equation [51 Thick solid lines indicate contours of 
constant J. 
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Fig. 12. — Upper limits on {av) as a function of calculated using Equation [TO] with a 
composite neutralino spectrum and the J values obtained for the Draco (solid line) and Ursa 
Minor (dashed line) NFW models and the M15 NFW+AC model (dotted line). Shown as 
open circles are MSSM models that fall with in ±3 standard dev iations of the relic density 
measured in the three-year WMAP data set (ISpergel et al.l 120071 ). 



